Flux dynamics and vortex phase diagram of the new superconductor MgB2 
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Magnetic critical current density and relaxation rate have been measured on MgB2 bulks from 
1.6 K to Tc at magnetic fields up to 8 Tesla. A vortex phase diagram is depicted based on these 
measurement. Two phase boundaries Hl^^^iT) and Hf^^(T) characterizing different irreversible flux 
motions are found. The H^^^^ (T) is characterized by the appearance of the linear resistivity and is 
attributed to quantum vortex melting induced by quantum fluctuation of vortices in the rather clean 
system. The second boundary Hf^^iT) reflects the irreversible flux motion in some local regions due 
to either very strong pinning or the surface barrier on the tiny grains. 

74.25. Bt, 74.20.Mn, 74.40.+k, 74.60. Ge 

The recently discovered new superconductor MgB2 
generates enormous interests in the field of superconduc- 
tivity Many important thermodynamic parameters 
have already been derived, such as the upper critical field 
7Jc2(0) = 13 - 20.4 T [|-|, the Ginzburg-Landau param- 
eter K « 26 1^, and the bulk critical superconducting 
current density jc « 2 x IQ^A/crn? at 4.2 K and 1 T ||] 
in bulk samples and jc ~ 8 x IQ'^A/crn? at 4.2 K and 
12 T in thin films. One big issue concerns however 
how fast the critical current will decay under a magnetic 
field and in which region on the field-temperature ( H- 
T ) phase diagram the superconductor can carry a large 
critical current density { jc )• This jc is controlled by 
the mobility of the magnetic vortices, and vanishes at 
the melting point between the vortex solid and liquid. 
A finite linear resistivity p = {E/j)j^Q will appear and 
the relaxation rate will reach 100% at this melting point 
showing the starting of the reversible flux motion. In 
this Letter we present an extensive investigation on the 
flux dynamics by magnetic relaxation method. A vortex 
phase diagram will be depicted based on the magnetic 
relaxation data. 

Samples investigated here were fabricated by both high 
pressure ( HP ) ( P = 6 GPa at 950° C for 0.5 hour ) 
and ambient pressure ( AP ) synthesis which was de- 
scribed very clearly in a recent publication ||^. High 
pressure synthesis is a good technique for producing the 
MgB2 superconductor since it can make the sample more 
dense and uniform ( in sub-micron scale ) and also pre- 
vent the oxidization of Mg element during the solid reac- 
tion. Our HP samples are very dense and look like metals 
with shiny surfaces after polishing. Scanning electron mi- 
croscopy ( SEM ) shows that the HP sample is uniform 
in sub-micron scale but some disordered fine structures 
are found in 10 nanometer scale, being similar to the in- 
ternal structure of large grains seen in the AP sample. 
X-ray diffraction ( XRD ) analysis on both type of sam- 
ples show that they are nearly in a single phase with the 
second phase ( probably MgO or MgBj^ ) less than 1 
wt.%. For simplicity we present in this Letter only the 
results from the HP sample. 



Resistive transition was measured by the standard 
four-probe technique and the magnetic measurements 
were carried out by a superconducting quantum interfer- 
ence device ( SQUID, Quantum Design MPMS 5.5 T ) 
and a vibrating sample magnetometer ( VSM 8T, Oxford 
3001 ). To precisely calculate the critical current density 
jc the sample has been cut with a diamond saw into a 
rectangular shape with sizes of 3.2 mm ( length ) x 2.7 
mm ( width ) x 0.4 mm ( thickness ). Fig.l shows the 
diamagnetic transition of one of the samples measured 
in the field-cooled ( FC ) and zero-field-cooled ( ZFC ) 
process. In the FC process, the temperature was lowered 
from above Tc to a desired temperature below Tc under a 
magnetic field, and the data are collected in warming up 
process with field. Its signal generally describes the sur- 
face shielding current and the internal frozen magnetic 
flux proflle. In the ZFC process, the temperature was 
lowered from above Tc to a desired temperature below 
Tc at a zero field and the data are collected in the warm- 
ing up process with a field. Its signal generally describes 
the internal magnetic flux profile which is ultimately re- 
lated to the flux motion. Resistive and diamagnetic mea- 
surement show that the transition is very sharp with a 
perfect diamagnetic signal and T""**^* = 40.3 K with the 
transition width only 0.4 K. 

In Fig.l we show the magnetization hysteresis loops ( 
MHL ) measured at temperatures ranging from 2 K to 38 
K. The symmetric MHLs observed at temperatures up to 
38 K indicate the dominance of the bulk current instead 
of the surface shielding current. The MHLs measured at 
low temperatures, such as 2 K, 4K show quite strong flux 
jump which will be discussed elsewhere A small kink 
appearing in the field ascending process in low field region 
may be attributed to the Bean-Livingston surface barrier 
and will be presented separately jl^] . It is easy to see that 
the MHLs measured at low temperatures ( e.g., 2 K to 
10 K ) are too close to be distinguishable. From these 
MHLs one can calculate jc via jc — 20AM/Fa(l — a/36) 
based on the Bean critical state model, where AM is the 
width of the MHL, V, a and b are the volume, length 
and width of the sample, respectively. The result of jc is 
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FIG. 1. Magnetization hysteresis loops measured at 2, 4, 
6, 8, 10, 14, 18, 22, 26, 30, 32, 35, 37, K and 38 K ( from 
outer to inner ). All curves here show a symmetric behavior 
indicating the importance of bulk current instead of surface 
shielding current. The MHLs measured at low temperatures ( 
e.g., 2 K to 10 K ) are too close to be distinguishable. Strong 
flux jump has been observed at 2 K and 4 K near the central 
peak. 

shown in Fig. 2. It is clear that the bulk critical current 
density jc of our sample is rather high. For example, at 
T = 18 K and ^iqH = 1 T, we have jc = 1.1 x lOM/cm^, 
which is among the highest values that have been re- 
ported P,pd| in bulk samples. The jc{H) curves have 
been measured with three different field sweeping rate 
200Oe/s, 100 Oe/s and 50 Oe/s. It is interesting to note 
that for all MHLs, there are small tails ( shown in the in- 
set to Fig. 2 ) in high field region. Accordingly the jc{H) 
curves show also a small tail in high field region. Since 
jc value in the small tail region is almost 4 x 10^ times 
smaller than that in the low field region, therefore it is 
safe to conclude that this small tail is due to some sec- 
ondary effect, such as some local regions with very strong 
pinning or the surface pinning by the tiny grains. From 
the contour of jc vs. H shown in Fig. 2 one can see that 
there are two ways to determine the so-called irreversibil- 
ity line Hirr- The first way is to take a criterion for jc 
— 30 A/cm'^ just before the appearance of the tail. This 
is named as the bulk irrever sibility fine H^rr^{T) which 
signals the sharp drop of jc before the tail and reflects 
the irreversible flux motion in clean bulk samples. The 
second way is just to follow the small tail to a criterion 
of ?>A/cm^ . In this way one can determine a higher irre- 
versibility line Hf^^ which marks a boundary above which 
a complete reversible flux motion occurs. As mentioned 
above, the small tail in our present samples is due to a 
secondary effect, therefore the major part of the vortex 
system melt at H^^^^. Thus we may conclude that the 
H^^^^ reflects the melting of the vortex matter in the 
clean limit, for example, in single crystals. If more pin- 
ning centers are introduced into the sample, one should 
be possible to push the IL more close to the upper critical 
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FIG. 2. Critical current density jc calculated based on the 
Bean critical state model. At each temperature the data has 
been measured with three field sweeping rate: 200 Oe/s, 100 
Oe/s, 50 Oe/s. The faster sweeping rate corresponds to a 
higher dissipation and thus higher current density. From these 
data one can calculate the dynamical magnetic relaxation rate 
Q. The jc{H) curves measured at low temperatures are very 
close to each other showing a rather stable value of Hf^''' 
when T approaches zero K. The dashed horizontal line marks 
the criterion of jc = SOA/cm^ for determining the J/f"/*. 
While by following the small tail ( which has been attributed 
to a secondary effect ) to a criterion of jc = 3A/cm^ , another 
phase line -fff^^ is determined. 

fleld Hc2{T). The " irreversibility lines " by following 
these two different methods are shown in Fig. 3. It is 
clear that the lower phase boundary H^^^^ terminates 
at about 8 T at 2 K, while the extrapolated values ( 
Hf^^ ) reaches about 12 T at 2 K. There is a large sep- 
aration between these two phase boundaries. It is im- 
portant to note that Larbalestier et al. |jl^ have found 
the similar feature where they regard the lower phase 
line H^^^^ as the Kramer line. Fig. 3 shows also the up- 
per critical fleld Hc2{T) determined from the tempera- 
ture dependent magnetization by deflning the Hc2{T) as 
the point at which the magnetization start to deviate 
from the normal state linear background |l^ and that 
determined from the resistive measurement by Takano 
et al. 1^ on HP samples. A striking result is that the 
H^^l^{T) extrapolates to a rather low field, here for ex- 
ample, H^^^^{Q) « 8 T, while the Hc2{T) extrapolates 
to a much higher value ( w 15 T ) at zero K. There 
is a large separation between the two fields 77c2(0) and 
Hirr{0). If following the hypothesis of the vortex liquid 
above Hirr, we would conclude that there is a large re- 
gion for the existence of a quantum vortex liquid at zero 
K. This can be attributed to a quantum fiuctuation effect 
of vortices in bulk MgB2- Although the lowest temper- 
ature in our present experiment is 1.6 K, however, from 
the experimental data one cannot find any tendency for 
Hirr'^iT) to turn upward rapidly to meet the Hc2{0) at 
zero K. One may argue that the Hirr{T) probably can 
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FIG. 3. H-T phase diagram for the new superconductor 
MgB2- The circles represent the bulk irreversibility lines 
Hi^''' : open circles, this work measured by VSM; filled 
circles, another sample measured by VSM; half-filled circles, 
HP sample measured by SQUID. The triangles represent the 
Hf^^iT) : filled triangles, determined directly from the small 
tail; open triangles, following the tendency of tail ( shown 
by the 3 dashed lines ) to a criterion of jc = ZA/crn^. 
The squares represent the upper critical field Hc2{T): filled 
squares, Takano's data from resistive measurement; open 
squares, this work from the M(T) measurement by SQUID. 
All solid lines are guide to the eye. 

be increased to higher values by introducing more pinning 
centers into this sample. This is basically correct since re- 
cently a higher irreversibility line Hirr{T) has been found 
in thin films 0| and bulk samples irradiated by heavy ions 
I p^ . This indicates that the quantum vortex melting in 
our present bulk sample is due to the strong quantum 
fluctuation which smears the perfect vortex lattice lead- 
ing to the vanishing of the shear module Cge of the vor- 
tex matter ( probably within grains ). Dense disorders 
will strengthen the shear module and thus enhance the 
bulk irreversibility line. The irradiation of heavy ions by 
Perkins et al. did not suppress but strongly increase 
jc at a high field would suggest that the low value of 
H^^^^{T) measured in unirradiated bulks ( like our sam- 
ples here ) is not due to the weak links since otherwise 
the jc value would drop even faster with increasing the 
magnetic field after the irradiation. This indicates that 
the rather low Hf^''J:'^{T) observed in our present work 
reflects the melting of the vortex matter in sample with 
relatively pure structural details rather than the breaking 
of the Josephson couplings between grains. In addition, 
as mentioned above, the H^^^'^iT) found in present work 
is very close ( or identical ) to those by other authors 
|pT|,p[ found in bulks, therefore it shows a more intrinsic 
feature corresponding to the melting of vortex matter in 
the rather pure system. All these imply a strongly linked 
current flow between grains ]l^ . 

In order to investigate the flux dynamics in the vortex 
solid state below H^^J:'' we have carried out the dynamical 



FIG. 4. Field dependence of the relaxation rate at temper- 
atures of 2, 4, 6, 8, 10, 14, 18, 22, 26, 30, 32, 35, 37 K. The 
dashed line is a guide to the eye for 2 K. It is clear that Q 
will rise to 100% at 8 T at 2 K. Since Hf"''' is rather stable at 
low temperatures, it is anticipated that Hi^^''{0) ~ 8T being 
much smaller than Hc2{0) « 15 T. 

relaxation measurement. According to Schnack et al. 
JlSf and Jirsa et al. jl^], in a held sweeping process, 
if the field sweeping rate is high enough, the quantity 
Q — dlnAM / dln{dH / dt) is identical to the relaxation 
rate S — —dlnM/dlnt measured in the conventional re- 
laxation method, where Q is called as the dynamical re- 
laxation rate, AM the width of the MHL, dH/dt the 
field sweeping rate. The raw data with three different 
sweeping rate ( 200 Oe/s, 100 Oe/s, and 50 Oe / s ) are 
shown in Fig. 2. The Q values vs. field for different tem- 
peratures are determined and shown in Fig. 4. It is clear 
that the relaxation rate increases monotonically with the 
external magnetic fleld and extrapolates to 100% at the 
bulk melting point H^^J:'' . A closer inspection to the data 
flnds that the rising of the relaxation rate Q vs. fleld is 
slowered down when it runs into the small tail region. 
Since the small tail is a secondary effect as mentioned 
before, we leave it to a future investigation. Here we 
concentrate on the flux dynamics before the setting in of 
the small tail. At 2 K it is found from the Q(II) data 
that the melting fleld is about 8 T, being very close to 
Hf;;"^-(T = 2K) determined from the jc(ff) curve. It 
is known that the H'j'^'-^^{T) is rather stable in low tem- 
perature region, therefore we can anticipate a rather low 
value of i/f"^"^(0) which is below 9 T being much lower 
than i/c2(0). The large separation between Hf^^''{0) and 
_ffc2(0) may manifest the existence of the quantum vortex 
liquid due to strong quantum fluctuation of vortices in 
the pure system of MgB2. Fig. 5 shows the temperature 
dependence of the relaxation rate Q. The arrows point 
at the irreversibility temperatures at the corresponding 
flelds H^^^^{T). It is clear that the relaxation rate in 
wide temperature region keeps rather stable against the 
thermal activation and fluctuation. However, when the 
bulk melting point H^rr^{T) is approached the relaxation 
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FIG. 5. Temperature dependence of the relaxation rate at 
fields of 1, 2, 3, 4, 5, 6 T. The dashed lines are guides to the 
eye. The arrows point at the bulk irreversibility temperatures 
T^^^^ determined from the curve. 

rate will quickly jump to 100%. This may indicate that 
the thermal fluctuation is not the dominant process for 
the flux depinning in the new superconductor MgB2. It 
shows high possibility for the vortex quantum melting 
even at a finite temperature. Worthy to note is that the 
quantum tunneling rate is extremely low at an intermedi- 
ate field, such as Q = 0.7% at 2 K and 2 T, but is rather 
high at a high field, for example Q — 33% at 2 K and 6 
T. This may imply that the field will greatly enhance the 
vortex quantum fluctuation and tunneling. 

In conclusion, in rather pure samples of MgB2, the 
irreversibility field is rather low comparing to the upper 
critical field in low temperature region showing the possi- 
ble existence of the quantum vortex liquid due to strong 
quantum fiuctuation. The temperature and field depen- 
dence of the relaxation rate may further suggest that the 
vortex melting at a finite temperature is also induced by 
strong quantum fluctuation in pure systems, such as sin- 
gle crystals and bulks. A small tail following the jc{H) 
curve in high field region is attributed to a secondary ef- 
fect induced by strong pinning or surface barrier on some 
tiny grains. This opens a way for the higher irreversibil- 
ity line when more pinning centers are introduced into 
the pure sample. 
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